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Abstract 

We analyze the Wilson loop for a pure Yang-Mills theory, using a decoupling solution in close 
agreement with lattice computations. At one-gluon exchange level it is seen that the potential 
cannot yield a linear rising contribution as expected for a confining theory. Next-to-leading order 
" correction gives rise to a quartic term for momenta in the gluon propagator that, in agreement 

q; with Gribov's view, yields a linear confining term. This correction is due to a two-loop or sunrise 

integral that we need to evaluate in the low-momenta limit. In the infrared regime, the physical 
consistency of the theory is determined by a natural cut-off, arising from the integration of the 
classical equations of the theory, fixing in this way the regularization scheme. 
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I. INTRODUCTION 

Since the groundbreaking work of Kenneth Wilson [lj , a proof of confinement for quantum 
chromodynamics in the continuum limit is lacking. This is expected to be one of the major 
goals for our understanding of low-energy behavior of hadronic matter. Since then, a lot of 
fundamental work has been pursued reaching a deeper comprehension of the problem 2|-[5]. 

Stated in its simplest way, the confinement problem can be expressed through the be- 
havior of the expectation value of the Wilson loop that, for a free theory, is just reduced 
to the computation of a double path integral on the propagator of the theory. When this 
expectation takes the simple form exp(— <r r S), being S the area spanned by the loop and a r 
a constant called string tension, one has that the area law holds and the potential rises lin- 
early with the distance granting completely bounded states. This result has been confirmed 
by computations on lattice {2], lj] and is generally stated as "Cornell potential" being the 
sum of Coulombian and linearly rising terms. 

The main problem to face to understand confinement is that the theory at low energies 
is strongly coupled. This means that perturbation theory is no more reliable to perform 
computations. Besides, very few techniques are available to manage the theory in this limit 
and, generally, one resorts to numerical studies on a lattice. Indeed, recent studies of the 
propagators of a pure Yang-Mills theory produced some striking results. These results on 
huge volumes, arriving to such a significant value as (27/m) 4 , were presented at the Lattice 
2007 Conference in Regensburg 6MS] • The shocking conclusion was that the scenario devised 



since then 



11] , generally accepted as correct, was not describing the situation seen on 



the lattice in three and four dimensions: The gluon propagator was reaching a plateau at 
lower momenta with a finite non-zero value in zero, rather than being zero there, and the 
ghost propagator was behaving as that of a free massless particle while an enhanced one 
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was expected. All in all, the running coupling was seen to bend clearly toward zero 
without evidence of a non-trivial fixed point as was expected instead. Anyhow, there is not 
a widespread agreement on the definition of the running coupling in the infrared limit. 

The scenario that emerged has an interesting interpretation. Assuming that an infrared 
trivial fixed point for a pure Yang-Mills theory exists 13Hl6|. the theory can be easily 
managed through a perturbation expansion. The propagator we obtain can be used to get 
several properties of the theory and QCD as well [171 ] . In this paper, we will show that this 
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propagator, by itself, implies a screened potential and so, it needs a next-to-leading order 
correction to recover a Cornell potential and confinement. It is important to point out that, 
while Yang-Mills theory has an infrared trivial fixed point, this is no more true for QCD 
due to the presence of quarks. Besides, what makes small a running coupling is just the 
decreasing of momenta that are expected to go to zero like the fourth power. As we will see, 
the excellent agreement with numerical data, till a perfect coincidence for numerical solution 
of Dyson-Schwinger equations [3 1 ls strongly supporting this view. This implies that 
the definition of the running coupling discussed in these papers is acceptable in the infrared 
limit and we are exploiting this possibility. 

Finally, we emphasize that a decoupling solution for the gluon propagator does not grant 
a potential behaving as the Cornell potential with a linearly rising term but rather a screened 
potential. This has been recently shown, numerically solving Dyson-Schwinger equations, by 



Gonzalez, Vento and Mathieu 



20] . As we will show in this paper, the reason for this is that 



we need loop corrections to the propagator in the infrared limit. In this way, a quartic term 
in the momenta is seen to appear, in agreement with a Gribov-like propagator, providing 
the needed linearly rising term in the potential. 

The paper is structured as follows. In secHIl we present the formalism applied to scalar 
field and Yang-Mills theories. In sec JIIIt we show a comparison with numerical data strongly 
supporting the formalism. In sec JIVj we give the next-to-leading order correction to the 
propagator. In secE] we discuss the renormalized gluon propagator we obtained. In sec I VII 
we compute the interquark potential and the string tension. Finally, in sec JVIIl we present 
our conclusions. 



II. INFRARED QUANTUM FIELD THEORY 



In the eighties, Carl Bender and others 2jJ, |22j proposed a new approach to cope with a 



strong coupled quantum field theory. This pioneering view was too radical and so, produced 
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16. 



23|, 



too singular results to be useful. But one can improve on it as we showed recently 
assuming a bare coupling formally going to infinity. As we will show, immediate conclusions 
can be obtained both for scalar and Yang-Mills quantum field theories. 
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A. Scalar field theory 



We consider a massless neutral scalar field with the following generating functional 



Z[j] = N J [d(f>) exp U J d A x 



\{d4>) 2 -- A 4> 4 + j<l> 



(i) 



In order to get a strong coupling expansion, we formally rescale space-time coordinates as 
x — > yXx. Then, the functional takes the form 



Z[j\ = N J [d<f>] exp jz J d 4 x 



2A V V! V A" 



(2) 



But j is an arbitrary function and so we rescale it as jj A — > j. Now, a perturbative expansion 
can be written down as 

oo 

= 5>~"0 n (3) 



n=0 

that yields the following terms into the action 

50 = [ d 4 x 

51 = 

So = / d 4 x 



i(90 o ) 2 -^+J0o 



+ d(j) d(f)2 



^2 + j<j>2 



(4) 
(5) 
(6) 



d 4 x [<90 o <90i - 0o0i + 3<f>i] 

>.)'-§< 

We note that this expansion requires that we are able to solve the classical equation of 
motion 

d 2 0o + A0 o i =j (7) 
in the limit A — > oo. Indeed, making use of this equation the following functional is obtained 



(8) 



Here 0o is a functional of the current ?. In this form, we are not able to manage this 



expression. So, we make the Ansatz 
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25) 



0o = A* / d x'A(x — x')j(x') + 



and we need to solve the following equation for the Green function 



d 2 A(x - x') + A[A(x - x')f = -5\x - x'). 



(9) 



(10) 



Here /i is an arbitrary constant with the dimension of energy arising in the solution of the 
classical theory. Then, using an iterative procedure, the next-to-leading order correction is 



/i\ J d 4 x'A(x-x') j/i J d 4 x"[A(x' -x")fj(x"^ 



d 4 x"A(x' — x")j(x") 



+ .... 



(11) 



Higher order corrections can be similarly obtained. Finally, the propagator can be obtained 



explicitly 



being 



23 



261 ] . One has 



B r 



n=0 



V 



m l n + ie 



Br, 



(2n + i; 



7T 



2 (_l)« e -(n+|h 



K 2 (i) l + e -(2n+l)7r ' 

with K(i) = dO/Vl + sin 2 6 « 1.3111028777 and a mass spectrum 



(12) 



(13) 



m. 



1 

n + - 



7T 



A\ 3 



(14) 



2 y V2. 

We recognize that the theory develops a mass gap due to the fmiteness of the coupling. 
This result would be very difficult to obtain with a weak coupling expansion. But, another 
striking result is that this theory is infrared trivial as the generating functional in this limit 
is just Gaussian 



Zo\j] = iVexp 



^ / d 4 xd 4 yj(x)A(x - y)j(y) 



(15) 



This can be also evinced from the form of the propagator (1121) that presents just simple poles 
and can be ascribed to a free theory of massive particles with a superimposed spectrum of 
a harmonic oscillator. Finally, the next to leading order term takes the form 



Z\j]nZ [j] /[#]exp<{- / d'x 



(16) 



B. Yang-Mills theory 



The foundations for a possible understanding of QCD at low-energies were posed on the 
eighties [27 



28] . People realized that, with a current expansion, a sensible low-energy limit 
was obtained. But not a great step forward was possible due to the stumbling block of an 
unknown gluon propagator. This was a major obstacle to complete this program. We will 



now see how a consistent scenarion can indeed be built on this basis. We have the generating 
functional 

Z[j,e,e] = iVexpj-z J d 4 x ^TrF 2 + ^{d ■ A) 2 + {(fd^c a + g^f abc d^c c ) jx 



exp 



J d 4 x (j;^ a + e° c ° + cv) 



;i7) 



being F£ v = d^A^—d v A a -\-gf abc A b A'l the field strength and A a ^ the vector potential. Exactly 
in the same way as for the scalar field, we operate a rescaling of space-time variable through 
the 't Hooft coupling y/Ng* for a SU(N) gauge group. So, again, we have to solve the 
classical equations of motion 



(18) 



for an arbitrary gauge fixed through the parameter £. Now, we would like to study quantum 
field theory for Yang-Mills equations assuming that a trivial infrared fixed point exists. In 
order to get such a theory, we have to show that a set of classical solutions exist to build 
up the corresponding quantum field theory. Indeed, this is the case and such solutions are 



instantons 



291 ] . This can be proved quite easily by putting A a = rj^cp^x) being rf some 



constant coefficients. Then, Yang-Mills equations collapse to 



1 



(19) 



In the Lorenz (Landau) gauge, this equation is further simplified and we can use the fact 
that rjlif - = N 2 — 1 and so we set j<f> = r]^j ua . In this case, the theory reduces to that of 
the scalar field and we can easily do quantum field theory for a trivial infrared fixed point. 
Quantum field theory will be build up on the instanton field that is replicated for each 
component of the Yang-Mills field. 

Now, we are able to give explicitly the gluon propagator in the Landau gauge as 



D*(p) = 5 ab (v, - 2^) A(p) + O (-±=- 



(20) 



with A(p) given by eq.f ll2j) . Here A — > Ng 2 and we interchange the constants \l and A. So, 
we can write the spectrum for SU(N) as 



m r 



1 

n+ 2 



71 



K(i) 



Ng' 



A. 



(21) 



This can be proven very easily from the definition of the two-point function. We have 

D^(x-y) = (TAl(x)A!t(y)). (22) 
Using our instanton solutions, one sees that 

K( x ~v) = vZvl(mx)(f>(y)) + o(i/VNg) 

= r ] ;r ] b u A(x-y) + 0(l/VNg) (23) 

giving our result. A possible choice for components of the gauge field is 

rU = S ab [rj, v - ^) (24) 

being r]^ the Minkowski metric. We are now able to develop a quantum field theory for the 
Yang-Mills field supporting a trivial infrared fixed point. 

For the sake of completeness, we give here the ghost propagator. Ghost field just decouples 
in the infrared limit due to the mapping theorem and so, 

G(p) = ^— + 0(l/^g). (25) 

Then, the generating functional for the gauge field in the infrared limit takes the following 
Gaussian form 



Z[j] = N'exp 



i J d 4 xd 4 yj^(x)D^(x-y)f b (y) 



+ ° [ttJ (26) 

entailing triviality. The spectrum has a mass gap with a superimposed spectrum of a har- 
monic oscillator as if such particles would have a structure. We note that infrared triviality 
does not apply to QCD due to the presence of quarks. 



III. COMPARISON WITH NUMERICAL DATA 



A priori, it is not clear why our classical solutions should be selected to build up a 
quantum field theory. Even if our approach is consistent, we cannot be sure that something 
should enter, and it is instead neglected, into the solutions needed for the theory. So, we 
need to compare our conclusions with numerical data from lattice and numerical Dyson- 
Schwinger equations. These studies have been pursued for some decades and reached their 
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maturity with a significant breakthrough on 2007 where the gluon propagator was seen to 
reach a finite non-zero value lowering momenta js-8]. This was also confirmed by numerical 
studies of Dyson- Schwinger equations [18 ]. 

Firstly, the most important point to be assessed is the existence of the trivial fixed point 
at infrared (momenta going to zero). As already said, this can be seen from recent lattice 



computations [12( and numerical solution of Dyson-Schwinger equations [18J, [19] provided 
one agrees on the definition of the running coupling at lower momenta. Assuming this, we 
will show below the behavior of the running coupling for the latter case (see fig. IIIII) . We 
note also that, given a different definition of running coupling, a similar result has been 
given in |3o| . 




FIG. 1: Running coupling obtained by numerically solving Dyson-Schwinger equations (here a 
-0.5059 ± 0.0012, b = 1.045 ± 0.001 and c = 0.004 ± 0.001. 



The definition of the running coupling is the one considered in 
It is given by 

a s (p 2 ) = aMZ{p 2 )H\p 2 ) 



121 ] and proposed in |31|. 



(27) 



having introduced the dressing functions for the gluon and ghost propagators respectively 
as Z(p 2 ) = p 2 A(p 2 ) and H(p 2 ) = p 2 G(p 2 ). In our fit from numerical solution to Dyson- 
Schwinger equations, it is seen to be very near the zero value as momenta goes to zero. 
This result is identical to the one presented in 12j. We note that, with this definition, the 



8 



running coupling goes to zero like p 2 while, with the definition given in |30j, it runs to zero 
like p 4 . 

Below we give the comparison with numerical results for the propagators. Data for 80 4 
and Dyson- Schwinger equations were given directly from measurement while measurements 
at 128 4 recovered from fig. 2 in [Tj. One sees that, increasing volumes, our propagator tends 
to the measured one and the agreement improves neatly in the deep infrared. Different mass 
gaps are due to different values of (3. For numerical Dyson-Schwinger equations [19!] , one has 
no volume problem being already in the infinite limit and, indeed, we get perfect coincidence 
of our propagator with the measured one. Note that we consider a weak dependence on the 
gauge group as showed in [32j that is fully consistent with our discussion above. 






p=6 (48) 4 


+ 


p=6 (64) 4 


* 


p=6 (80) 4 




P=6.2 (48) 4 


o 


P=6.2 (64) 4 


□ 


p=6.2 (80) 4 




- Frasca 




p 2 [GeV 2 ] 



FIG. 2: Gluon propagator in the Landau gauge for SU(3), 80 with a mass gap of mo = 321 MeV 
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FIG. 3: Gluon propagator in the Landau gauge for SU(2), 128 4 with a mass gap of nig = 555 MeV 
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FIG. 4: Gluon and ghost propagators in the Landau gauge for SU(3) obtained by numerically 
solving Dyson-Schwinger equations and a mass gap mo = 399 MeV 
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These numerical results give a strong support to our scenario and to the conclusion of 
the existence of a trivial infrared fixed point for Yang-Mills theory. 



IV. TWO-LOOP CORRECTION 

In sec. Ill B\ we have shown that, for a Yang-Mills theory with a trivial infrared fixed 
point, it is enough to study the corresponding quantum theory of the instanton solutions 
represented by a massless scalar field with a quartic self-interaction. For our aims, we 
need just to evaluate the next-to- leading order correction to the propagator and derive the 



potential from it as shown in |20j. This will move us away from the infrared fixed point and 
will provide a comprehension on the theory approaching such a fixed point. 
So, from eq.( !TT|) we have to evaluate 

A R (x - x') - A(x - x') =/j 2 J d 4 x"A(x - x") [A(x" - x')} 3 (28) 

where A(x — x') is the one given in eq. (fT2|) . We point out here that, for a Yang-Mills theory, 
we have to set A = Ng 2 . Turning to momentum space this gives 

A^> - W = AlA« / ^L^L E -J^ x (29) 

E B n2 sr^ B n3 



„ 2 Pi - m l 2 (P - Pi - P2) 2 - m 2 n3 



This integral is well-known in quantum field theory, it arises as a two loop sunrise diagram 



331 ] and, in the form given above, has not a closed form value. But, we are in the infrared 
limit and we need its value just in the limit of small momenta. But before we take small 
momenta limit, consistently with our approximation technique, we need to take the limit 
A — > 00. This will give immediately 

d^pi d 4 p 2 



A R (p 2 ) - A(p 2 ) = A-^A(p 2 ) B n,B n2 B n , j 



x 



y-v" ^ -m-wnaj ( 27r )4( 2vr V- 
' ni, 712,^3 \ / \ / 

1 1 

A§ (2n~i + l) 2 (2n 2 + l) 2 (2n 3 + l) 2 mg 

\_ ( (P-Pi- P2) 2 

A 2 \{2 ni + l) 2 (2n 2 + l) 2 (2n 3 + l) A ml 

+ pi , p\ 



{2m + l) 4 (2n 2 + l) 2 (2n 3 + l) 2 mg (2m + l) 2 (2n 2 + l) 4 (2n 3 + l) 2 m; 
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0/ J 

+ 0(A-*), (30) 
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where we have set, using eq. ([Tlj) . m = jjn / K(i) ps \l. These integrals are really singu- 
lar but we note that we cannot extend the integration range to infinity as we expect our 
computation to hold just for momenta p < fi, the cut-off of the theory that emerged by 
integration of the classical equations. So, the theory has a natural cut-off. In this way we 
can immediately evaluate the correction with the following integrals: 

d 4 p 1 d 4 p 2 



1 



(2tt) 4 (2tt) 



d A pi d 4 p 2 
(2tt) 4 (2tt) 4 



7T 8 



(p - pi - p 2 ) = -r p + -/i 



G? 4 P1 C? 4 p 2 2 

(2vr) 4 (2vr) 4Pl ' 2 



7T 
3 7T 8 



(31) 



So, we finally get 



A R (p 2 ) - A(p 2 ) = A(p 2 



1 27 13.3-48 



A3 7T 6 



A 7T 8 



1 



16 fl 2 



O A-S 



(32) 



We note that the cut-off disappears everywhere except as a scale for the momenta contri- 
bution as it should. Now, we can interpret this as an expansion in inverse of a/A of the 
denominator of the propagator and this means that we have to evaluate 



P 



1 27 13.3-48 
+ 

At 7T 



1 " 8 ' A 7T 8 



+ 



(P 2 

13.3-9p 4 



! 1 27 13.3-48 

\\ 7T 8 A 7T 8 

1 27 13.3-48 



16 /i 2 



A2 vr 



+ 



1 - 8 ' A 7T 8 



1 + 



3_p^ 
16 /x 2 



(33) 



A 7T 8 /i 2 

We get a mass renormalization factor causing the appearance of a p 4 term in agreement with 
Gribov expectations and, finally, looking at the mass term we see that 



1 27 1 3.3-48 

AItt 8 



A 



7T 



1 + 1^ 

16 fi 2 



(2n+l) 2 m 2 X^ 



( 1 27 13.3-48 

\k 7T 8 A 7T 8 



1 + — ^ 



3_p^ 

16 /x 2 
(34) 



and we get a renormalized running coupling 



K(p 2 



1 - 



1 27 13.3-48 



AU 8 



3_^_ 
16 /i 2 



A 7T 8 

From this we get the renormalization constant for the field 



O \~i 



(35) 



Up 2 ) 



1 27 13.3-48 



A^ vr 8 



+ 



A vr 8 



1 + 



16 /i 2 



O A"* 



(36) 
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V. RENORMALIZED GLUON PROPAGATOR 



Once we obtained the field renormalization constant Z<p(p 2 ), we can provide a next-to- 
leading order formula for the gluon propagator. Let us introduce the constant Z = Z^(0) 
and A = Ng 2 . We can rewrite the gluon propagator in the Landau gauge in the following 
form 

On 



where we have set 



ab I T]fiu 



m 2 n (p 2 



P 



Z + 



t , 1 3.3-9 p 4 \ ^2 ( n 2\ 

13.3-9p 2 



CM A - * 



O A"2 



(37) 



(38) 



A 7T 8 /i 2 

We note that the coefficients B n are exponentially damped with n and so, just the first few 
terms give the physics of infrared Yang-Mills theory. This means that our pr opagato r agrees 



fairly well with the one obtained in the Refined Gribov-Zwanziger scenario [34 



. So, we 



can attempt to evaluate the condensates after the correction is introduced. As already done 
in [37)], we consider just two terms in our propagator and write down 



G(p 2 



Z 1 B 



Z Q x Bx 



In order to estimate the condensates, we note that Bq + B\ w 1 obtained from the exact 
relation B n = 1 accounting for the exponential damping of higher coefficients. Now, by 



+ 



1 3.3-9 p 4 



+ m\ (p 2 



(39) 



the definition given in [38 

(g 2 A 2 ) 



of the condensate, we get 
9 iV 2 — 1 



13 N 



Z [(m 2 + mi) - (B mj + B x ml)\. 



(40) 



In our previous work 



37| we get about 0.13 GeV 2 for this condensate and our result here is 



not so much different from it. Again, this is a proo 



our scenario agrees fairly well with that given in 



of existence of such a condensate and 



3J-|36|. 



VI. WILSON LOOP AND POTENTIAL 



In order to compute the potential in a pure Yang-Mills theory at the infrared fixed point, 
we have to evaluate 

^^^^ = /Wc^ 



J[dA][dc\[d 



c e 4 



J d 4 xTrF 2 +iS g [c 
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being S g [c, c] the contribution of the ghost field, T a the anti-hermitian generators of the 
gauge group and V the path ordering operator. In our case, in the infrared limit, we have 
a trivial fixed point. This implies that our generating functional takes a Gaussian form 



Z [j] = exp 



d A x'd i x"j a ^x')D a hx' - x")j bv (x") 



(42) 



and so, the evaluation of the Wilson loop is straightforwardly obtained as [4 



W[C] = exp 



~C 2 (R) j> dx» j dy u D^(x - y) 



(43) 



being C 2 (R) the quadratic Casimir operator that for SU(N) in the fundamental represen- 
tation, R — F, is 02(F) = (N 2 — 1)/2N. In our case, the propagator has the form given 
in secJTT] for the Landau gauge. The fall-off to large distances of this propagator grants 
that ordinary arguments to evaluate the integrals on the path apply. Indeed, using Fourier 
transform one has 

.2 , d 4 p 



W[C] = exp 



~C 2 (R) 



(27T) 



-A(p ) V/j,v - 



PlJ.Pi 
p2 



J j> dx" j> dy u e- ip( - x - y ) 



We need to evaluate 



1(C) = rtfa, 



j) dx^ j) 



dy u e 



v -ip(x-y) 



(44) 



(45) 



provided the contributions coming from taking into account the term run faster to zero 
at large distances. This must be so also in view of the gauge invariance of Wilson loop. 
With the choice of time component in the loop going to infinity while distance is kept finite, 
we can evaluate the above integral in the form 



1(C) w 2vrT5(p )e" 



ipx 



(46) 



and we are left with 



This yields 



being 



W[C]^e W -T 9 ^C 2 (R) J ^A(p,0)e-^ 



W[C] = exp [-TV YM (R)\ 



V YM (R) = - 9 ^C 2 (R) j -0-A(p,O)e-**. 



(47) 



(48) 



(49) 
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This equation should change if we consider a running coupling but, for the moment, we will 
not do this as does not change too much our conclusions. So, in our case we have to evaluate 



the integral 20] 



2 poo 

V YM (R) = -^jf*( R ) J o dp P A(p 2 )sm(pR) (50) 

using eq.( l37|) . If we neglect the contribution of the sunrise diagram we see that the propaga- 
tor is a sum of Yukawa propagators. This means that the above integral will yield a sum of 
Yukawa potentials and we get an overall screened potential that does not grant confinement 
as discussed in [20j. So, we see that in the deep infrared, one-gluon exchange does not grant 
a confining theory. But, accounting also for our two-loop correction, we will have to compute 
instead 

Vym (R) = -^C 2 {R)Z, 1 Y B n r dp P^(P R ) . (51) 

8nR ° ^ r.2 I 1 3.3-9 p 4 , m 2/^ V ; 

We can now evaluate this integral by rewriting it in the form 

v ; 8nR K ' 3.3-9^ y (p 2 + K 2 ) 2 -K 4 K ' 

where we have considered the limit A — > oo, vr? w p 2 = 0(y/\) even after a renormalization 
of mass — > Z^vr? n is taken, B„ = 1 and we have set k 2 = Kw^-- This integral can be 



computed exactly to give 



2 

Vym(R) « -^C 2 {R)e~^ R sinh(-|it!) (53) 



that gives at small values of kR 

2 T 7T 7T _ 

(54) 



Vym(R) ~ ~/^C 2 (R) 
8irR 



^=kR--k 2 R 2 + 0((kR) 3 ) 
y/2 2 U ' ' 



So, remembering that A = d(r)g 2 , being d(R) = N for SU(N) in the fundamental represen- 



tation, we get a linear rising potential with a = jj-C^k as expected. We can rewrite this 



4 4tt 

as 



7T-2 2 / C 2 (R)d(R) f . 



This result should be compared to the d — 2 + 1 case [41] that has 



, jC 2 (R)d(R) 



V^a+i ~ 9 \l ^ • (56) 



The agreement is strikingly good. 
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VII. CONCLUSIONS 



We have shown that, in the framework of our formalism, confinement is achieved for QCD 
satisfying the area law from the expectation value of the Wilson loop once next-to-leading 
order correction is accounted for. We started from the hypothesis of the existence of a trivial 
infrared fixed point for a Yang-Mills theory and this drove us to the fundamental conclusion 
that the scenario generally depicted by Wilson and Gribov is essentially correct. 

As an aside, it interesting to note as our approach recovers in some limit the refined 
Gribov- Zwanziger scenario confirming the existence of a A 2 condensate. This condensate is 
relevant for the existence of a mass gap in the theory. 

The overall agreement with lattice data and numerical studies of Dyson-Schwinger equa- 
tions gives a serious hint that this view can represent an important track toward a full 
understanding of quantum field theory in the infrared limit. 
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